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Finite Size Effects in Quark-Gluon Plasma Formation 
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Using lattice simulations of quenched QCD we estimate the finite size effects present when a gluon plasma 
equilibrates in a slab geometry, i.e., finite width but large transverse dimensions. Significant differences are 
observed in the free energy density for the slab when compared with bulk behavior. A small shift in the critical 
temperature is also seen. The free energy required to liberate heavy quarks relative to bulk is measured using 
Polyakov loops; the additional free energy required is on the order of 30 — 40 MeV at 2 — 3 T c . 



Lattice gauge theory has given us considerable 
information about bulk thermodynamics, but fi- 
nite size effects have up to now been studied using 
simplified, phenomenological models. We study 
the behavior of a gluon plasma restricted to a 
slab geometry, with the longitudinal width much 
smaller than the transverse directions. This in- 
ner region is heated to temperatures above the 
bulk deconfincmcnt temperature, surrounded by 
an outer region which is kept at a temperature 
below the deconfinement temperature, providing 
confining boundary conditions for the inner re- 
gion. Details of this work are given in Q . 

By allowing the coupling constant (3 = 6/g 2 
to vary, a spatially dependent temperature can 
be introduced. We have chosen the temperature 
interface to be sharp, in such a way that the lat- 
tice is divided into two spatial regions, one hot- 
ter and one colder. The quenched approximation 
simplifies the role of the cold region, because be- 
low T c , the dominant excitation at low energies 
is the scalar glueball. Temperatures near T c are 
smaller than glueball masses by about a factor 
of four, so glueballs play no the pressure in the 
hadronic phase is essentially zero. We thus expect 
the slab thermodynamics to be largely insensitive 
to the precise temperature of the region outside 
the slab, as long as it is sufficiently low. In full 
QCD, this insensitivity to the outer temperature 
would not hold, due to pions. Note that the role 
of boundary conditions here is quite different from 



that in bubble nucleation. In that case, both (3i n 
and p ou t are taken to be near T c . [|||| 

The slab is given a fixed lattice width, w = 6a, 
rather than of fixed physical width. Since N t = 4, 
wT is fixed at 3/2. At higher temperatures, the 
width in physical units is somewhat smaller than 
the longitudinal size of the plasma formation re- 
gion expected in heavy ion collisions. While the 
use of equilibrium statistical mechanics to study 
gluon plasma properties during the early stages of 
plasma formation may appear suspect, a simple 
estimate using the Bjorken model [Q shows that 
when a coin-shaped region of width 1 fermi has 
expanded to 1.5 fermi, the variation in tempera- 
ture is only from 0.8TJ) at the center of the coin 
to To at its edges. 

The free energy density / for the slab was ob- 
tained using the standard method || of integrat- 
ing the lattice action with respect to j3. We use 
a convenient convention for the sign of / that is 
opposite the usual one. In the bulk case, / is then 
identical to the pressure p. 



J_\0 

t a 1/3. 







dp [(S)t - (S) ] 



Ha 



(1) 



*We thank the U.S. Department of Energy for financial 
support under grant number DE-FG02-91-ER40628. 



where S = (l/N)ReTr U p . As in the bulk case, 
it is necessary to subtract the zero-temperature 
expectation value from the finite temperature ex- 
pectation value, in this case using the same pair 
of j3 values. In general, quantum field theories 
with boundaries develop divergences that are not 
present in infinite volume or with periodic bound- 
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ary conditions. It is expected that any countert- 
erms would be gauge-invariant and have dimen- 
sion < 3; there are no such counterterms. || 

In order to take advantage of the data on bulk 
thermodynamics provided by the Bielefeld group 
, we worked consistently with lattices of overall 
size 16 3 x 4. The values used for each subtraction 
come from 16 4 lattices with identical values of /3j„ 
and /3 ut- The value of f3 ou t was held fixed at 5.6 
while Pi n varied from 5.6 to 6.3. For comparison, 
the bulk transition for N t = 4 occurs at P c (N t = 
4, N s = 16) = 5.6908(2) p c (N t = 4, N s = oo) = 
5.6925(5). § 

Figure 1 shows the free energy density f/T 4 
versus T/T c compared with the bulk pressure. 
The free energy in the slab is lower than the bulk 
value by almost a factor of two at 2T C . It appears 
that the slab value is slowly approaching the bulk 
value, but other behaviors are also possible. Cal- 
culations of the finite-temperature contribution 
to the Casimir effect for a free Bose field con- 
tained between two plates show that f/T 4 has a 
non-trivial dependence on the dimensionless com- 
bination u>T.|| H It is natural to ask if the cor- 
rections to the free energy seen here can be ac- 
counted for by the conventional Casimir effect. A 
straigtforward calculation of the free energy of a 
non-interacting gluon gas confined to a slab shows 
an increase in the free energy density over the 
bulk value by a factor of about 1.63 at wT = 3/2. 
The Casimir effect alone is thus unable to explain 
the reduction of the free energy observed in our 
simulations. We are currently investigating more 
elaborate theoretical models which include the ef- 
fect of a non-trivial Polyakov loop. 

A consistency check was performed on the sur- 
face effects. JToJ The free energy density was cal- 
culated for a system at Pi n = 6.0 by perform- 
ing simulations with /3j„ fixed at 6.0 and ft out 
varying from 5.6 to 6.0. Combining these re- 
sults with the bulk data of reference |?J creates 
a path equivalent to varying (3i n while holding 
l3 ou t fixed. For (3 in =6.0 and (3 ou t =5.6, this 
gives f/T 4 = 0.65 ± 0.04, to be compared with 
f/T 4 = 0.69 ± 0.03 for the direct calculation. 

The major source of systematic error lies with 
the choice of boundary conditions for the slab, 
here set by Pout- We have estimated the system- 
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Figure 1. Slab free energy compared with bulk 
pressure. 



atic error associated with the choice of bound- 
ary conditions by by performing simulations at 
Pin = 6.2 and p ou t = 5.5 on 16 3 x 4 and 16 4 
lattices. These results suggest that lowering p ou t 
from 5.6 to 5.5 reduces the free energy by roughly 
10 percent at Pi n = 6.2. 

We define an effective surface tension a(w,T) 
by 



/ = p — 2a(w,T) /w 



(2) 



where the notation a(w,T) recognizes that the 
surface tension a does depend on the width of the 
slab and the internal and external temperatures. 
In the limits where T approaches T c and w goes 
to infinity, this quantity approaches a^. Figure 
2 shows a{w,T)/T 3 versus T/T c for wT = 3/2; 
representative error bars are shown. Away from 
the bulk critical point, a/T 3 rises quickly to a 
peak at about 1.4T C , and then falls slowly as T 
increases. A large non-equilibrium surface ten- 
sion has also been observed in measurements of 
the equilibrium surface tension ao, where these 



effects were obstacles to obtaining qp. 10 

The multiplicative divergence in the Polyakov 
loop can be eliminated when comparing bulk ex- 
pectation values to those in finite geometries. We 
define 



AF Q {x) = —T In [PslabW/Pbulk] 



(3) 
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Figure 2. Non-equilibrium surface tension versus 
T. 



as the excess free energy required to liberate a 
heavy quark in the slab geometry relative to bulk 
quark matter at the same temperature. 

In figure 3, we show the expectation value for 
the Polyakov loop versus z measured in lattice 
units for several values of (3. Each curve is nor- 
malized by dividing the values of the Polyakov 
loop by the bulk expectation value at the corre- 
sponding value of p. Error bars are shown only 
for even values of z. It is clear that a significant 
change occurs between (3 = 5.8 (T = 1.23 T c ) 
and = 5.85 (T = 1.36 T c ). For larger values 
of f3, AFq diminishes to a value of approximately 
30 - 40 MeV in the middle of the slab. 

There are significant deviations in the slab ge- 
ometry from bulk behavior and the ideal gas law, 
arising from a strong non-equilibrium surface- 
tension. This non-equilibrium surface tension can 
be an order of magnitude greater than the equilib- 
rium value. Surface tension effects also produce 
a mild elevation of the apparent critical tempera- 
ture. Measurement of Polyakov loop expectation 
values relative to bulk shows that the suppres- 
sion of heavy quark production due to the slab 
geometry is small. 

There are good reasons to call our result an 
estimate rather than a calculation. Although lat- 
tice gauge theory simulations of bulk behavior can 
be made arbitrarily accurate in principle, in this 



case there is uncertainty in the exterior bound- 
ary conditions appropriate; the applicability of 
equilibrium thermodynamics at this early stage of 
quark-gluon plasma formation can also be ques- 
tioned. However, this result is the best estimate 
available now, and further refinements are possi- 
ble. 
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Figure 3. Polyakov loop versus z. 
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